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The origin of Gamma-Ray Bursts is one of the most interesting puzzles in recent 
astronomy. During the last decade a consensus formed that long GRBs (LGRBs) arise 
from the collapse of massive stars and that short GRBs (SGRBs) have a different origin, 
most likely neutron star mergers. A key ingredient of the Collapsar model that explains 
how the collapse of massive stars produces a GRB is the emergence of a relativistic jet 
that penetrates the stellar envelope. The condition that the emerging jet penetrates the 
envelope poses strong constraints on the system. Using these constraints we show that: (i) 
Low luminosity GRBs (//GRBs), a sub population of GRBs with a very low luminosities 
(and other peculiar properties: single peaked, smooth and soft) cannot be formed by 
Collapsars. //GRBs must have a different origin (most likely a shock breakout), (ii) On 
the other hand regular LGRBs must be formed by Collapsars. (iii) While for BATSE 
the dividing duration between Collapsars and non-Collapsar is indeed at ~ 2 sec, the 
dividing duration is different for other GRBs detectors. In particular most Swift bursts 
longer than 0.8 sec are of a Collapsar origin. This last results requires a revision of many 
conclusions concerning the origin of Swift SGRBs which were based on the commonly 
used 2 sec limit. 
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1. Introduction 

Gamma ray bursts (GRBs) are among the most amazing transients known. In a 
few second a GRB emits the energy that a star like our sun emits in its whole 
life time. Their origin has puzzled astronomers since their serendipitous discovery 
in the late sixties. After two decades in which it was believed that the GRBs are 
Galactic it was realized , in the early n i neties, that they have a cos mological origin 
[Meegan et al., (1992)[ |Piran (1992)] |Mao fc Paczynski (1992)] . The distance 



scale set immediately the energy scale to be ^ 10 5i erg (including beaming 
corrections that were realized t owards the end of the nineties |Rhoads (1997) [ 
Rhoads (1999)| |Sari et al.(1999) | ) . Together with the short time scale this has led 



inevitably to the conclusion that the events involve the formation of a newborn 
compact object, most likely a black hole. This conclusion has left practically just 
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two progenitor candidates: A collapsing massive star or the merger of two neutron 
stars (or a neutron star and a black hole). 

The observations of a few (long) GRB afterglows in 1997 revealed that those 
bursts arose in star forming regions. Paczynski (1998) who noticed that, quickly 
suggested that long GRBs (LGRBs) are related to core collapse events. At roughly 
the same time MacFadyen & Woosley (1999) suggested the Collapsar model. 
According to this model a GRB is produced by a relativistic jet that emerges 
from the center of a massive collapsing star and penetrates the stellar envelope. 
By now the name Collapsar is used with different variation^] We stress that here 
we use this original definition for a Collapsar: A jet that penetrates the envelope 
of a Collapsing star. Using numerical simulations MacFadyen & Woosley (1999) 
demonstrated that a relativistic jet can indeed penetrate a stellar envelope. Again 
roughly at the same time Galama et al. (1998) discovered that GRB 980425 
was associated with the powerful type Ic supernova: SN 1998bw. However GRB 
940425 was a strange GRB. It was very weak, with an energy a few orders of 
magnitude less than the energy of a typical GRB. Additionally it was single peaked 
and smooth and it had a very soft spectrum. It was not clear that this was a 
regular GRB and hence the association of GRB 980425 with SN 1998bw was 
not sufficient to demonstrate a GRB-SNe association. Shortly after that Bloom 
et al. (1999) and others discovered red bumps in the afterglows of more distant 
aregular GRBs. These red bumps were interpreted as the signatures of 1998bw- 
like SNe, supporting the GRB-SNe association. However the evidence for a GRB- 
SNe association was inconclusive until SN 2003dh was discovered in association 



with the regular GRB 030329 |Hjorth et al., (2003) | |Stanek et al., (2003)] . Since 
then a few other GRB-SNe associations were discovered. Even though most of 
these GRB-SNe associations are with weak, smooth, single peaked GRBa^]this is 
generally considered as a "proof" of the Collapsar model for LGRBs. 

An inspection of B ATSE's GRBs' temporal distribution revealed 
|Kouveliotou et al., (1993)"] two group s: short (Tqq < 2 sec) and long (Tqq > 2 sec.). 
Already in 1995 it was pointed out |Cohen fc Piran (1995) Piran (1996)] that 



the two groups have a different spatial distribution. The observed short GRBs 
(SGRBs) are significantly nearer (and weaker). This suggested the possibility of 
different physical origin for the two populations. As it takes time (and energy) 
to cross the relatively large stellar envelope it was argued that SGRBs cannot 
be produced by Collapsars |Matzner (2003)] . In most cases Collapsars produce 
LGRBs, but by now we know that in some cases Collapsars produce SGRBs (see 
£4 below). However a variant on this original argument that we discuss here (in 
£3 and El below) shows that most SGRBs cannot be produced by Collapsars. 
Lack of detection of SGRB afterglows left the situation inconclusive until 2005, 
when Swift localized the first short bursts and the first SGRBs' afterglows were 
detected. It turned out that SGRBs are not associated with star forming regions 
(some arise in elliptical galaxies) and as such they are not associated with deaths 
of massive stars. The progenitors could be neutron star mergers (as suggested 

In some cases the term Collapsar is used generically for any model that involve a collapsing 
star, regardless whether there is a jet or not. In other cases it is used more restrictively to imply 
a situation in which the collapsing star produces an accreting black hole as a central engine that 
drives a relativistic jet. 

2 Another association of a regular regular GRB and an SN, GRB 101219B and SN 2010ma, 
was discovered recently. 
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already in 1989 |Eichler et al., (1989)|). I 
demonstration of this origin |Nakar (2007)] 



owever as yet there is no conclusive 



We describe here new results, derived by Bromberg et al. (2011a, 2011b, 
2012a, 2012b) concerning the nature of GRB progenitors. We briefly discuss, in 
^2] some recent analytic results concerning relativistic jet penetration through the 
stellar envelope [Bromberg et al., (2011a)l . We then consider their implications on 



this picture. In c|3]we demonstrate that GRBs, those that appear in most GRB - 
SNe associations, cannot be produced by Collapsars Bromberg et al., (2011b)|. 
While this weakens the case for the association of regular LGRBs with SNe, 
we show in £j4] that when combined with the GRBs' temporal distribution 
these considerations dem onstrate that the LGRBs originate from Collapsars 
|Bromberg et al., (2012a)] , providing a direct observational indication for jets that 
puncture the stellar envelope. Further inspection of the temporal distribution 
enables us to estimate (in q5J), for the first time, t he fraction of Collapsars am ong 
SGRBs as a function of the observed duration jBromberg et al., (2012b)] . We 
show that this fraction depends strongly on the detector (in particular on its 
spectral window). In particular the standard limit of 2 sec is invalid for Swiff s 
observations, for which a limit of 0.8 sec is much more appropriate. 



2. Jet Propagation 

A schematic picture of a relativistic jet propagating within a stellar envelope is 
depicted in Fig. [T] There are a few critical com ponents. A a double shock system 



appears at the head of the jet Matzner (2003)] . While the jet is highly relativistic 
these shocks slow down the head and it typically propagates with a sub or mildly 
relativistic velocity. The hot material that streams sideways out of the jet's head 
produces a cocoon that engulfs the jet. While expanding sideways into the rest 
of the stellar envelope (this expansion will eventually blow out a fraction of the 
stellar envelope) it a lso squeezes the jet and prod uces a (radiative) collimation 
shock within the jet | Bromberg fc Levinson (2007)] . 

As long as the jet is within the stellar atmosphere all its energy is dissipated 
at the jet's head. The total dissipated energy equals therefore the jet's luminosity 
times the time it takes to cross the envelope. Since the inner engine is much 
smaller than the envelope it is decoupled from the jet that crosses the envelope 
on a much larger scale and one can expect that the luminosity before and after the 
jet breaks out are comparable. Using the observed GRB luminosity to estimate 
the jet p ower before breakout we can estimate the duration of the dissipation 
phase as |Bromberg et al., (2011a)] : 

/ £. \-l/3 / Q \ 2/3 / ^ \ 2/3 / ^ \ 1/3 

t B *15 -c^ io51 j;° /sec j ^— j (^-j , (2.1) 

where Lj so is the isotropic equivalent jet luminosity, 6 is the jet half opening angle 
and we have used typical values for a LGRB. R* and M* are the radius and the 
mass of the progenitor star, where we have normalized their value according to the 
typical radius and mass inferred from observations of the few SNe associated with 
LGRBs. For the jet to break out, the central engine must continue operating for 
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a duration longer than If the inner engine stops before the jet's head crosses 
the envelope the jet won't produce a regular GRB. 



3. Low luminosity GRBs 

The duration of the prompt emission, approximated by T90, is given simply by: 

Tgo = t e — tB, (3.1) 

where t e is the total time that the engine powering the jet is active. Within the 
Collapsar model, without fine tuning, only a small fractions of the bursts should 
have Tgo/tB <C 1 (see Q. Namely, it is unlikely that the engine operates just long 
enough to let the jet break out of the star and then stops right after breakout. This 
argument was used by Matzner (2003) to argue that Collapsars cannot produce 
SGRBs, for which Tgo/iB^l- This is indeed confirmed in Fig [2] in which the 
distribution of Tgo/is is shown for both LGRBs and SGRBs. One can clearly see 
two distinct populations LGRBs for the majority of which Tgo/ts > 1 and SGRBs, 
all of which satisfy Tgo/ts < 1. 
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Figure 2. The distribution o f T 90 /t B for LGRBs , //GRBs and SGRBs (from 

[Bromberg et al., (2011b)] ). 



Fig [2] depicts also a third group of GRBs, low luminosity (//GRBs). Like 
SGRBs, the observed duration distribution of //GRBs is inconsistent with the 
predictions of the collapsar model. In particular a large fraction of //GRBs have 
Tgo/ts <C 1. The probability that the observed //GRBs Tgo/is distribution is 
consistent with the LGRBs distribution is smaller than 5% implying that //GRBs 
cannot be generated by Collapsars and they must have a different origin. 

//GRBs are a group of six GRBs whose luminosities are around 10 47 - 10 49 
ergs/sec, at least two orders of magnitude below the average luminosity of a typical 
GRB. Remarkably //GRBs are not characterized just by their low luminosity. 
They are also single peaked, smooth and soft. //GRBs include GRB9890425 
(the first GRB detected to accompany a Supernova - 1998bw) as well as 
a few other GRB-SN pairs: GRB 031203/SN20031w; GRB060218/SN2006aj; 
GRB100316D/SN2010bh. GRB051109B shows all the common features of //GRBs 
but it lacks a repor ted SN. It is associat ed with a star forming region in a spiral 
galaxy at z = 0.08 [Perley et al., (2006)] . All //GRBs are at very low redshifts. 
With such low luminosities they couldn't have been detected from further out. 
While only six //GRBs have been observed so far, given these distances, the //GRBs 
inferred rate per unit vol ume is much larger than the rate of regular LGRBs 
|Soderberg et al., (2006a)] . In fact this rate is so large that //GRBs cannot be 
significantly beamed as even with a modest beaming corrections the rate would 
have exceed the rates of their associated SNe - broad line type Ibc. 
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An interesting and likely possibility is that //GRBs' jets are weak and 
fail to break out from their progenitors. A"failed jet" dissipates all its energy 
into the surrounding cocoon and drives its expansion. As the cocoon reaches 
the edge of the star its forward shock may become mildly or even ultra 
relativistic emitting the 7-rays observed in //GRBs when it breaks out. 
This idea that //GRBs arise from shock breakouts was suggested shortly 



following the observations of GRB980425/SN1998bw [K ulkarni et al., (1998). 
MacFadyen, Woosley, fc Heger (2001)] |Tan, Matzner McKee (2001)| . It drew 



much more attention following the observation of additional //GRBs with similar 
properties and especially with t he observation of a th e rmal component in 
the spectrum of Z/GRB 060218 |Campana et al., (2006) [ |Wang et al., (2007) [ 



|Waxman, Meszaros, fc Campana (2007)| . Yet, it was hard to explain how shock 
breakout releases enough energy in the form of 7-rays. Katz, Budnik & 
Waxman (2010) realized that the deviation of the breakout radiation from 
thermal equilibrium provides a natural explanation to the observed 7-rays. More 
recently, Nakar & Sari (2012) calculated the emission from mildly and ultra- 
relativistic shock breakouts, including the post breakout dynamics and gas- 
radiation coupling. They find that the total energy, spectral peak and duration of 
all //GRBs can be well explained by relativistic shock breakouts. Moreover, they 
find that such breakouts must satisfy a specific relation between the observed 
total energy, spectral peak and duration, and that all observed //GRBs satisfy 
this relation. These results lend a strong support to the idea that //GRBs are 
relativistic shock breakouts. From a historical point of view this understanding 
closes the loop with Colgate's (1968) original idea, that preceded the detection of 
GRBs, that a SN shock breakout will produce a GRB. 

As we discuss in the following section, the observed GRB duration distribution 
indicates the existence of many "failed jets" in which the engine time is shorter 
than the breakout time. This is consistent with the observations that the rate of 
//GRBs is much higher than the rate of regular LGRBs. 



4. Long GRBs and Collapsars 

As most of the GRBs associated with SNe are //GRBs one might think at first that 
this new understanding rules out the Collapsar model for LGRBs. However, on 
the contrary, these arguments provide a new and unexpected direct observational 
confirmation of the Collapsar origin of LGRBs. Consider again Eq. |3.1| Under 
very general conditions this equation results in a flat duration distribution for 
durations significantly shorter than the typical breakout time. 



It follows from Eq. 3.1 that the distribution, p 7 (Tgo) of the observed GRB 
durations is a convolution of p e (t e ), the distribution of engine operating times, 
and pb^b) the distribution of jet breakout times. Under quite general conditions 
(more specifically unless p e varies very rapidly around ts, an unlikely situation) 
the following limits hold: 



^ (T90) ~\ Pe(T 90 ) T 90 »t B ■ 
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Particularly interesting for our purpose here is the appearance at short 
durations X90 <J,i_B of a flat region, a plateau, in which the rate of events is 
independent of the duration. Remarkably such plateaus exist in all the observed 
GRB duration distributions (see Fig. [3]). They weren't no ticed so far because 
the "canonical" distribution plot KouveTiotou et al., (1993)] depicts dN / dlog(T) 
instead of dN/dT. 

The higher end of the plateau enables us to estimate ts and from this to 
infer some the basic properties of the collapsing stars. We find, for example, 
that a typical progenitor size is ~ 5R@. Another interesting feature seen in Fig. 
[3] is th e ra pid decline at durations longerr than ts- At this regime according 
to Eq. - 



4.1 



the distribution is dominated by p e (t e ), thus p e (Tgo) ~p 7 (Tgo). An 
extrapolation of this distribution to shorter engine operating times suggests that 
there are numerous cases in which t e <ts and the jet fails to break out. This is 
in a very nice agreement with the very large inferred event rate of //GRB if these 
are interpreted as"failed jets". 
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Figure 3. The duration distributions, dN/dT go , of BATSE (red), Swift (blue) and Fermi GBM 
(green) GRBs. Also plotted is the distribution of the soft (hardness ratio < 2.6) BATSE bursts 
(magenta). For clarity the Swift values are divided by a factor of 5 and the Fermi GBM by 
15. Note that the quantity dN/dT is depicted and not dN / dlogT as traditionally shown in 



such plots |Kouveliotou et al., (1993) . The black lines show the best fitted flat interval in 
each data set: 5-25 sec (BATSE), 0.7- 21 sec (Swift), and 2.5 - 31 sec (Fermi). The upper 
limits of this range indicate a typical breakout time of a few dozens seconds, in agreement 
with the prediction of the Collapsar model. Soft BATSE bursts show a considerably longer 
plateau (0.4 — 25 sec), indicating that most of the soft short bursts are in fact Collapsars (from 
|Bromberg et al., (2012a)] .) 
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The appearance of these plateau (and their dependence on the observed 
hardness) is the first direct observational confirmation of the Collapsar model. 
More specifically these plateau demonstrates the fact that the duration of an 
LGRB is the difference between two independent time scales. This confirms a basic 
prediction of the Collapsar model: the overall duration is the difference between 
the time that the engine operates and the time it takes the jet to penetrate the 
stellar envelope. 

At very short durations the plateau doesn't extend all the way to zero. This 
does not rule out the model. At this regime non-Collapsar SGRBs, that have 
a different origin, appear and dominate the distribution. As different detectors 
have different relative sensitivities to long (and soft) vs. short (and hard) GRBs 
the duration at which short non-Collapsars begin to dominate varies from one 
detector to anothe^ To demonstrate this dependence on the detector's spectral 
window we artificially change BATSE's effectiveness for detection of hard SGRBs 
by considering only softer BATSE bursts (hardness ratio < 2.6). As expected, for 
this softer BATSE sample the non-Collapsar peak shrinks and the plateau extends 
down to shorter durations. 



5. Short non - Collapsar GRBs 

Shortly after Kouveliotou et al. (1993) demonstrated that there are two 
populations of GRBs: long and short ones it became clear that these 
two populations have different spa tial distributions and a different origin 
[Cohen fc Piran (1995)1 l piran (1996)"] . By now we know that LGRBs arise from 



Collapsars. SGRBs have other projenitors, most likely neutron star mergers. As 
we are still uncertain concerning the origin of SGRBs we denote them here as 
non-Collapsars. So far it was implicitly assumed that the division line between 
long and SGRB is at 2 sec regardless of the observing satellite. The existence of a 
plateau in the observed LGRBs (of Collapsar origin) duration distribution enables 
us to determine, for the first time, the fractions of Collapsars vs. non-Collapsars 
as a function of the observed time for every specific detector. While we cannot 
determine if a specific bursts is a Collapsar or not we can give now a probabilistic 
estimate for a given duration and hardness. 

The basic idea is very simple. For a given detector we determine the rate 
of detection of bursts within the plateau. This provides an estimate for the 
detection rate of short duration Collapsars by this detector. Now we can compare 
this rate to the rate of SGRBs at any given duration and obtain the Collapsar and 
non-Collapsar fractions as a function of duration. This estimate is performed for 
different detectors or even for different detection windows (hardness) for a specific 
detector. 

We have fitted the different duration distributions with a plateau (representing 
Collapsars) and a lognormal distribution (for non-Collapsars). The fit is 



3 In principle we should have worked with a redshift corrected sample. However such samples 
of SGRBs are too small and we are forced to use the observed durations. This imply that we have 
to worry also about time dilation which, in turn depend on the sensitivity of the detector and 
the corresponding depth of the samples. However, as SGRBs are detected only from relatively 
small redshifts this time dilation correction is not significant. 
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remarkably good and it enables us to estimate the fraction, f^c of non-Collapsars 
from the total number of observed GRBs as a function of the observed duration, 
T90 (see Fig. W|. For Batse T90 < 2 sec is a reasonable threshold to identify non- 
Collapsars. This limit results in a probability > 70% for a correct classification 
for BATSE bursts. However, this condition is misleading for Swift bursts. At 
T90 = 2 sec a Swift burst has a 84±14% probability to be a Collapsar! Clearly, for 
Swift a 2 sec division line results in a large number of misidentified Collapsars as 
non-Collapsars. We propose to draw the division line between collapsars and non- 
Collapsars at the duration where the probability that a GRB is a non-Collapsar is 
50%. With this condition, a BATSE GRB can be classified as a non-Collapsar if 
its T90 < 3.1±0.5 sec. Swift bursts can be identified as non-collapsars only if their 
duration T90 < 0.8±0.3 sec, while the corresponding limit for GBM is T90 < 1.71q5 
sec. The results shown here can be expanded and improved when we consider the 
hardness of the bursts. 

Using these results one can go back and examine various studies of short 
bursts that attempted to compare various features of short (standing for non- 
Collapsars) with long (Collapsars) and check the samples used. A preliminary 
inspection of such studies reveals that while in some cases the short sample has 
only a small number of potential Collapsars in other cases the short sample 
used was heavily contaminated by high probability potential Collapsars (with 
observed duration shorter than 2 sec) and it is possible and even likely that 
these have dominated the results. Note in particular that our results are not 
inconsistent with those of Berger (2011), who finds a larger dispersion in properties 
of Swift SGRBs (with T20 < 2 sec host galaxies and the positions of the bursts 
within the hosts, as compared with the more homogenous properties of LGRBs 
host and their position within the hosts. A mixed sample that contains both 
Collapsars and non-Colla psars is expected to show such larger dispersion (see 
jBromberg et al., (2012b)] for further details). 



6. Conclusions 

To conclude we summarize our basic findings. We define^] a Collapsar as a 
collapsing massive star that produces in it center a relativistic jet. The jet 
penetrates the stellar envelope and produces the GRB ones it has left the star. 
The duration of the envelope penetration phase depends on the jet's luminosity, 
its opening angle and the size and density profile of the stellar envelope. 

One can expect that the duration of a burst is typically comparable to or longer 
than the jet breakout time. Indeed, a comparison of the estimated jet breakout 
time of a typical LGRB shows that it is shorter than the observed duration of the 
burst. On the other hand the jet breakout time is much longer than the duration 
of a short burst. This provides the first indication that SGRBs are not produced 
by Collapsars. We have shown that a third group of low luminosity GRBs also 
don't satisfy this condition. This implies that //GRBs don't arise from Collapsars. 

Within the Collapsar model the observed duration of a GRB is the difference 
between the time its central engine operates, producing the jet, and the jet's 
breakout time. This directly implies that at short durations the rate of bursts 



4 Note that the term Collapsars has also been used in both wider and narrower contents. 
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Figure 4. From top to bottom: BATSE, Swift and Fermi GBM fractions, /jvc of non-Collapsars 
from the total number of observed GRBs as a function of the observed duration, Tgo. The shaded 
regions represent 67% confidence limits of /jvc- Also plotted in red are the Tgo values for which 
/ N c = 0.5 (from |Bromberg et al., (2012b)] ). 



produced by Collapsars should be independent of their duration. We have shown 
that such a behavior is observed in the duration distributions of all GRB satellites: 
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BATSE, Swift and GBM. This provides the a direct observational confirmation 
of the basic prediction of the Collapsar model and as such demonstrates the 
Collapsar origin of LGRBs. 

This last feature also enables us to determine the fraction of Collapsars within 
the observed SGRBs. This fraction depends on the characteristics of the detector. 
For BATSE the standard division between Collapsar and non-Collapsars is indeed 
at ~ 2 sec. However for the softer Swift many bursts shorter than 2 sec are of 
Collapsar origin. This might have led to some confusion in the past in interpreting 
observations of these short bursts as indications for properties of non-Collapsar 
GRBs. 
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